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The relation between the slope of the nuclear symmetry energy at sat- 
uration density and the neutron skin thickness is investigated. Constraints 
on the slope of the symmetry energy are deduced from the neutron skin 
data obtained in experiments with antiprotonic atoms. Two types of neu- 
tron skin are distinguished: the "surface" and the "bulk". A combination 
of both types forms neutron skin in most of nuclei. A prescription to cal- 
culate neutron skin thickness and the slope of symmetry energy parameter 
L from the parity violating asymmetry measured in the PREX experiment 
is proposed. 

PACS numbers: 21.10.Gv, 21.30.Fe, 21.65.Ef 

1. Introduction 

The symmetry energy is a quantity of a great importance in nuclear 
physics and astrophysics pQ El El HI [5]. ^ g° verns numerous isospin-dependent 
properties of nuclei such as the binding energy, the location of the drip lines, 
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the density distributions, as well as the reactions: giant resonances, heavy 
ion collisions, isospin diffusion, and multifragmentation. At the same time, 
the nuclear symmetry energy is crucial in the astrophysical calculations of 
neutron stars, supernova explosions and stellar nucleosynthesis. 

Despite of its importance in nuclear physics, the symmetry energy is 
hardly known at densities different from saturation [6l [7] . Nuclear models 
predict a large variety of its density dependence. The values of the im- 
portant parameter L, which describes the slope of the symmetry energy at 
saturation, are spanned over a wide range. The so called "soft" symme- 
try energy can be found for the D1S parametrization where L = 22 MeV 
whereas L = 118 MeV is calculated for the NL3 parameter set which pre- 
dicts " stiff' symmetry energy. To fix the symmetry energy at subsaturation 
densities is needed for a better description of the isospin dependence of 
nuclear interactions. 

The neutron skin thickness is one of the observables where symmetry 
energy shows up in the ground state of nuclei [8]. Neutron skin defined 
through the rms radii of protons and neutrons depends on the properties of 
the nuclear surface. The relative differences of the neutron and the proton 
distributions in this region are sensitive to the symmetry energy at the 
subsaturation densities. Thus important information about nuclear forces 
can be deduced from the analysis of the neutron skin properties. 

In this paper we will study how neutron skin depends on the slope of 
the symmetry energy at saturation by means of the droplet model (DM) of 
nuclei. The predictions of the parameter L of symmetry energy will be made 
based on the neutron skin measurements in the antiprotonic atoms. The 
non-negligible influence of the differences of the diffuseness of the neutron 
and the proton density distribution at the nuclear surface on the value of 
neutron skin thickness will be shown. The relations between the neutron 
skin thickness of 208 Pb and the slope of symmetry energy parameter L with 
the parity violating asymmetry measured in the PREX experiment [9] will 
be investigated. 

2. Determination of the slope of symmetry energy 

The energy per particle of asymmetric nuclear matter can be expressed 

as: 

e(p, 6) = e(p, 0) + c sym (p)5 2 + 0(5 4 ) , (1) 

where 5 = (p n — p p ) / p is the asymmetry between the densities of protons 
and neutrons. This formula defines the symmetry energy of a nuclear EOS 
c sym ( / o), which can be expanded around saturation in the form: 

Cs y m(p) = J-Le + \K syra e 2 + 0{e 3 ) , (2) 
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Fig. 1. (Color online) The neutron skin thickness in 208 Pb calculated in various 
models plotted as a function of the ratio J/ Q of the symmetry energy at saturation 
to the surface stiffness coefficient (left) and the slope of symmetry energy coefficient 
L (right). 



where e = (po — p)/(3po)- While the symmetry energy at saturation J = 
Csym(po) is constrained by the empirical information to be around 32 MeV 
[IQ1 [HJ E2 03] j the predictions for the slope L = 3pdc sym (p)/dp\ Po and the 
curvature K sym = 9p 2 d 2 c sym (p)/dp 2 \ P0 parameters are varying substantially 
among the different theoretical models [14] . 

The neutron skin thickness defined as the difference of the neutron and 
the proton rms radii 

Ar^ = <r 2 >y 2 - (rX /2 (3) 
is described in the droplet model [15] by the formula 
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The second term in this expression is due to the Coulomb repulsion of 
protons and the last term depends on the surface diffuseness of neutrons b n 
and protons b p . If as in the original papers of the droplet model b n and b p 
are assumed to be equal this term vanishes. The leading term in Eq. (4) is 
almost a linear function of relative neutron excess / = (N — Z)/A: 

3 J I-Ic 

2 ro Q i + $J/4Q)A-W ' 1 j 

where ^ = 2O0V3 • 

The droplet model shows a linear correlation of neutron skin in any 

heavy nucleus with the ratio J/Q of the symmetry energy at saturation 

to the surface stiffness coefficient. Such correlation for the case of 208 Pb 

calculated in several nuclear forces is plotted in the left panel of Fig. 1. 

More surprising is that the neutron skin in 208 Pb is also linearly correlated 

with the parameter L what has been observed previously [U [16] (see the 

right panel of Fig. 1). 

To explain this relation, first, we transform the expression (5) into the 

form 

t = fy [J ~ a sym (A)] A 1 / 3 (I - I c ) , (6) 

where a sym (A) = J /(l + (9J/4Q)A-V 3 ) is the symmetry energy coefficient 
of the finite nucleus in the droplet model. We have shown that in the heavy 
nuclei with mass A > 200, o sym (^4) is equal to c sym (p) at densities around 
0.1 fm~ 3 [17J. This relation is fulfilled in good approximation also for lighter 
isotopes if c sym (p) is taken at a slightly smaller subsaturation density pa [E] . 
Within the mean-field approach, such density can be parametrized in terms 
of the mass number of the finite nucleus as follows: pa = Po~ Po/ i^+cA 1 ^), 
with the constant c fitted to reproduce the density of 0.1 fm -3 for 208 Pb 
[see Ref. [T7] for further details]. In this way, a sym (A) in Eq. (6) can be 
replaced by c syra (pA) and using Eq. (2) we can obtain 

In the last formula the leading term is t = L I. As mentioned before, the 
coefficient J = c sym (po) of the symmetry energy at saturation density has 
been empirically constrained to be around 32 MeV [TOj, HH [T2], [13] , which is 
well reproduced by the successful nuclear forces [2]. However, the values of 
the parameter L of the slope of symmetry energy at saturation are varying 
more substantially among the different theoretical models [HI [TT] dU [19] . 
As J varies in a much narrower range than L in the nuclear models, and as 
if J is larger (smaller) also L is usually larger (smaller), Eq. (7) suggests 
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to leading order a linear correlation between the neutron skin thickness in 
208 Pb and the parameter L, that is confirmed by the theoretical calculations 
as it can be seen in Fig. 1. 

From the independent calculations of b n and b p in the semi-infinite nu- 
clear matter [144 fTTj the influence of the diffuseness term in Eq. (4) on the 
total neutron skin thickness can be estimated. In the case of 208 Pb we have 
obtained values from around 0.035 fm to around 0.055 fm in various nuclear 
models. These numbers mean that diffuseness term can contribute as much 
as 20% to 40% to the total neutron skin thickness and they are too large to 
be omitted. Moreover, the neutron skin thickness calculated in the droplet 
model including the diffuseness term is very close to extended Thomas-Fermi 
predictions, whereas without this term the two models differ significantly. 
These results show that the assumption of equal surface diffuseness of neu- 
trons and protons is not consistent with mean-field predictions of density 
distribution in a nucleus. 

The experimental neutron skins deduced from the antiprotonic atoms 
[20^ I2T] presented in Fig. 2 show a linear increase of neutron skin with 
/. This trend is consistent with the droplet model prediction in Eq. (4). 
Since we have shown that the slope of the neutron skin with I depends 
on the parameter L (see Eq. (7)), the value of this parameter may be 
constrained by experimental data. The least squares fit of the formula (4) 
to the experimental data from the antiprotonic atoms, also presented in 
Fig. 2, reproduces average increasing trend of measurements. From the 
fitting procedure including the diffuseness term, we have found the value of 
L = 55 db 25 MeV for the slope of symmetry energy, which points towards 
a " soft" nuclear symmetry energy. The uncertainty of the value of L comes 
from large experimental errorbars of the measurements. 

3. The "bulk" and the "surface" types of neutron skin 

Two basic scenarios leading to create neutrons skin can be distinguished. 
They are schematically shown in Fig. 3. In the first type, called "skin" 
mode, the neutron surface is shifted outwards the proton surface. The al- 
ternative scenario is an increasing of the diffuseness of the neutron surface 
relative to the proton one. It is called the "halo" mode. When two param- 
eter Fermi (2pF) function 



is assumed for the neutron and the proton density distributions, both types 
of neutron skin are easy to classify. The "skin" mode corresponds to 2pF 
parametrization with a n = a p , whereas the pure "halo" mode is obtained 
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I = (N-Z)/A 

Fig. 2. (Color online) Neutron skin measurements in antiprotonic atoms[3Ul [2T] as 
a function of neutron excess i\ Average line of the experimental data is also given 
as well as results of the fit of formula (4) to the average line. 




Fig. 3. (Color online) Schematic view of the "bulk" ("skin") and the "surface' 
(" halo" ) types of neutron skin. 



with C n = C p . The central density is not properly reproduced by the 
sharp-edge sphere of a radius C with the volume conservation condition 
|22j . Therefore, the aforementioned "skin" and "halo" definitions do not 
correspond to the droplet model expression for the neutron skin thickness 
and cannot be directly related to the particular terms in the neutron skin 
formula in Eq. (4). Instead, we propose to use the "bulk" and the "surface" 
contributions |2"2"] . They are based on the same idea as the " skin" and 
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Fig. 4. (Color online) Neutron skin calculated with the SkM* Skyrme force in the 
antiprotonic atoms (left) and in the Sn isotopic chain (right). The "bulk" and the 
"surface" contributions are also shown. 



the "halo" types but with special attention paid on preservation of the 
properties of the interior of a nucleus by the "bulk" part. This contribution 
is denned through the equivalent sharp-edge radii [23] of neutrons R n and 
protons R p 
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As it is clear from Eqs. (9) and (10), the "bulk" and the "surface" contri- 
butions can be also easily expressed by the 2pF parameters. 

In Fig. 4 the neutron skin thickness calculated with the SkM* parameter 
set is plotted for the nuclei investigated in the antiprotonic atom measure- 
ments as a function of neutron excess /. A linear increasing trend of neutron 
skin thickness can be easily noticed. The "bulk" and the "surface" contri- 
butions calculated for the same nuclei are also plotted in Fig. 4. To obtain 
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these values the mean-field density distributions have been parametrized by 
the 2pF profile. Special care was taken on the precise reproduction of the 
surface region. The "bulk" contribution roughly shows a linear increasing 
trend as the total neutron skin with a similar slope. In the light nuclei with 
small / this contribution is negative due to the dominant influence of the 
Coulomb repulsion. In the heavy nuclei it grows up and states roughly one 
third of the total value of neutron skin. The "surface" contribution is very 
scattered especially in the light nuclei. In the region of the large neutron 
excess this part remains at the same level and does not change significantly 
with I. It seems that single-particle effects in particular nucleus have an im- 
portant influence on the surface region of some nuclei what impacts strongly 
on the "surface" part of neutron skin. On the other side, the "bulk" contri- 
bution is governed by the classical liquid drop properties of nuclei and the 
influence of quantal effects is limited. 

In the Sn isotopic chain the single-particle influence on the nuclear den- 
sity distribution is limited to shell structure of neutrons. The increasing 
trend of neutron skin as well as both "bulk" and "surface" contributions 
can be noticed. Again neutron skin thickness is composed of the both con- 
tributions. Shell effects are pronounced at magic numbers (1 = and 0.24). 
A huge increase of the "surface" part can be noticed in the mid-shell nuclei 
whereas there is a kink with a local minimum for doubly magic nuclei. The 
"bulk" contribution shows up opposite and less pronounced behaviour. The 
"surface" part is larger than the "bulk" one in most of the nuclei of the tin 
chain. 



The PREX (Pb Radius Experiment) performed at Jefferson Lab aims to 
constraint the neutron radius of 208 Pb by means of parity violating electron 
scattering (PVES) [U [231 113 [26]. In this experiment polarized electrons 
are scattered on the 208 Pb target. Electrons couple with protons mainly 
via the electromagnetic (7) interaction and with neutrons via the weak 
interaction (Z°). For ultra-relativistic electrons, as in PREX conditions, the 
weak potential has opposite sign depending on the polarization (helicity) of 
the electron beam. Therefore, the elastic differential crosssection of right- 
handed (da^/dVt) electrons is different to the one measured for left-handed 
(d(T-/d£l) electrons. The difference of these cross sections allows one to 
define the parity violating asymmetry: 



4. Predictions of PREX experiment 
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Fig. 5. (Color online) Relation between A pv and Ar„ p calculated in various models. 
Predicted errorbars with arbitrary chosen central value is plotted with thin green 
errorbar and preliminary measured value [5] with thick red errorbar. 



It is measured at PREX at the fixed energy of 1.06 GeV and the scattering 
angle of about 5° (gi a b ~ 0.47 fm -1 ). We calculate the elastic differential 
cross sections of 208 Pb for a polarized electron beam by means of the dis- 
torted wave Born approximation (DWBA) [19]. For that we have modified 
the code used in Ref. [27J. In short, the DWBA consist on the exact phase 
shift analysis of the Dirac equation, 

[op + /3m e + V(r)]tl> e = E^ e (12) 

where ip e is the electron wave function and V(r) = Vc(r) =t Vw(r) is the 
electroweak potential felt by the ultra-relativistic electron beam depending 
on their helicity. The connection between the parity violating asymmetry 
and the neutron skin thickness of a nucleus can be qualitatively understood 
from the analytical expression of the former within the plane wave Born 
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approximation (PWBA) 



PWBA 



G F q 2 - 



4 sin 2 W + 



Fn(q) ~ F p (q) 
F p {q) 



(13) 



where the neutron and proton formfactors F n (q) and F p (q), respectively, are 
also written in PWBA. Therefore, in the low-momentum transfer regime, 

F(q) « 1 - q 2 (r 2 )/6 and A™ BA explicitly depend on the factor ({r 2 )\ /2 - 

{r 2 )]/ 2 ) x ((r 2 )l/ 2 + (r 2 )p~^ 2 ). For realistic results, however, one cannot use 
the simple PWBA and full DWB A calculations must be performed [9J [21] • 
In Fig. 5 the parity violating asymmetry in 208 Pb calculated in DWBA 
for a large sample of successful nuclear forces is plotted versus the neutron 
skin thickness calculated for the same parametrizations. Models that fail to 
predict the well-known charge radius and binding energy of 208 Pb are not 
allowed in Fig. 5. A very good linear correspondence of the parity violating 
asymmetry and the neutron skin thickness is found. From the measurement 
of the parity violating asymmetry, the neutron skin can be easily deduced 
using linear dependence of Fig. 5: 



From this relation, taking into account that all models properly repro- 
duce proton rms radius, one can notice that the expected 3% errorbars of 
A pv should allow to obtain the neutron radius of 208 Pb with 1% uncertainty. 
It is illustrated in Fig. 5 with the green errobars plotted at the arbitrary 
chosen value of A pv . The PREX experiment completed a successful run in 
2010. The statistical error of the first run of the experiment was larger than 



This preliminary value of A pv with the reached 10% accuracy, marked in Fig. 
5 by the thick red errorbar, does not indicate on the forces which properly 
describe both the parity violating asymmetry as well as neutron skin. The 
PREX-II proposal aims to improve the accuracy of the experiment [9]. 

The slope of symmetry energy parameter L can be also calculated for 
each force presented in Fig. 5. As A pv is linear with Ar np and Ar np is 
linear with L (see Fig. 1) we may also expect a linear relation between A pv 
and L. From the confrontation of these two quantities we can find a nice 
correlation with average fitted line |19j 



^(ppm) = 0.788 - 0.375Ar np (fm) 



(14) 



3% 



AT? = 0.6571 ± 0.0604(stat) ± O.OUO(syst) 



(15) 



A p „(ppm) = 0.750 - 0.000557L(MeV) 



(16) 



The PREX experiment can provide information not only about the neutron 
skin in 208 Pb but also the independent estimation of L. 
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5. Conclusions 

We have discussed various aspects of neutron skin and its relations to 
the symmetry energy of nuclear matter. The linear relation between the 
thickness of neutron skin and the slope of the symmetry energy coefficient 
L has been proven. With the use of this correlation we have deduced from 
the neutron skin measurements in the antiprotonic atoms relatively small 
value of L = 55 MeV which hints the soft character of symmetry energy. 

Neutron skin can be produced as a result of the shift of the neutron 
surface in relation to the proton one (the "bulk" type), or due to the dif- 
ferences in the surface thickness of proton and neutron matter in a nucleus 
(the " surface" type) . In the general case both mechanisms contribute to the 
creation of neutron skin. 

The precise prescription for calculating the neutron skin thickness in 
208 Pb and the parameter L form the parity violating asymmetry in the 
PREX experiment has been given. High linear correlation between A pv and 
A r np in Pb has been noticed. 
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